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ABSTRACT 

The existence of compact gravitational lenses, with masses around 0.5 Mq, has been 
reported in the halo of the Milky Way. The nature of these dark lenses is as yet 
obscure, particularly because these objects have masses well above the threshold for 
nuclear fusion. In this work, we show that they find a natural explanation as being the 
evolutionary product of the metastable false vacuum domains (the so-called strange 
quark nuggets) formed in a first order cosmic quark-hadron transition. 
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One of the abiding mysteries in the so-called standard 
cosmological model is the nature of the dark matter. It is 
nniversally accepted that there is an abundance of matter 
in the universe which is non-luminous, due to its very weak 
interaction, if at all, with the other forms of matter, ex¬ 
cepting of course the gravitational attraction. The present 
consensus (for a review, see (Turner 1999; Turner 2000)) 
based on recent experimental data is that the universe is 
flat and that a sizable amount of the dark matter is ’’cold”, 
i.e. nonrelativistic, at the time of decoupling (we are not 
addressing the issue of dark energy in this work). Specula¬ 
tions as to the nature of dark matter are numerous, often 
bordering on the exotic, and searches for such exotic matter 
is a very active field of astroparticle physics. In recent years, 
there has been experimental evidence (Alcock et al. 1993; 
Aubourg et al. 1993) for at least one form of dark matter 
- the Massive Astrophysical Compact Halo Objects (MA¬ 
CHO) - detected through gravitational microlensing effects 
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proposed by Paczynski (1986) some years ago. To date, there 
is no clear picture as to what these objects are made of. In 
this work, we show that they find a natural explanation as 
leftover relics from the putative first order cosmic quark - 
hadron phase transition. 

Since the first discovery of MACHO only a few years 
ago, a lot of effort has been spent in studying them. Based 
on about 13 - 17 Milky Way halo MACHOs detected in the 
direction of LMC - the Large Magellanic Cloud (we are not 
considering the events found toward the galactic bulge), the 
MACHOs are expected to be in the mass range (0.15-0.95) 
Mq, with the most probable mass being in the vicinity of 
0.5 Mq (Sutherland 1999; Alcock et al. 2000), substantially 
higher than the fusion threshold of 0.08 Mq. The MACHO 
collaboration suggests that the lenses are in the galactic 
halo. Assuming that they are subject to the limit on the 
total baryon number imposed by the Big Bang Nucleosyn¬ 
thesis (BBN), there have been suggestions that they could 
be white dwarfs (Fields et al. 1998; Freese et al. 2000). It is 
difficult to reconcile this with the absence of sufficient ac¬ 
tive progenitors of appropriate masses in the galactic halo. 
Moreover, recent studies have shown that these objects are 
unlikely to be white dwarfs, even if they were as faint as 
blue dwarfs, since this will violate some of the very well 
known results of BBN (Freese et al. 2000). There have also 
been suggestions (Schramm 1998; Jedamzik 1998; Jedamzik 
& Niemeyer 1999) that they could be primordial black holes 
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(PBHs) ( ~ 1 Mq ), arising from horizon scale fluctuations 
triggered by pre-existing density fluctuations during the cos¬ 
mic quark-hadron phase transition. The problem with this 
suggestion is that the density contrast necessary for the for¬ 
mation of PBH is much larger than the pre-existing density 
contrast obtained from the common inflationary scenarios. 
The enhancement contributed by the QCD phase transition 
is not large enough for this purpose. As a result a hue tuning 
of the initial density contrast becomes essential which may 
still not be good enough to produce cosmologically relevant 
amount of PBH (Schmid et al. 1999). Alternately, Evans, 
Gyuk & Turner (1998) suggested that some of the lenses are 
stars in the Milky Way disk which lie along the line of sight 
to the LMC. Gyuk & Gates (1999) examined a thick disk 
model, which would lower the lens mass estimate. Aubourg 
et al. (1999) suggested that the events could arise from self- 
lensing of the LMG. Zaritsky & Lin (1997) have argued that 
the lenses are probably the evidence of a tidal tail arising 
from the interaction of LMC and the Milky Way or even a 
LMC-SMC (Small Magellanic Cloud) interaction. These ex¬ 
planations are primarily motivated by the difficulty of recon¬ 
ciling the existence of MACHOs with the known populations 
of low mass stars in the galactic disks. 

Adopting the viewpoint that the lensing MACHOs are 
indeed in the Milky Way halo, we propose that they have 
evolved out of the quark nuggets which could have been 
formed in a first order cosmic quark - hadron phase transi¬ 
tion, at a temperature of ~ 100 MeV during the microsec¬ 
ond era of the early universe. The order of the deconfine¬ 
ment phase transition is an unsettled issue till date. Lattice 
gauge theory suggests that in a pure {i.e. only gluons) SU (3) 
gauge theory, the deconfinement phase transition is of hrst 
order. In the presence of dynamical quarks on the lattice, 
there is no unambiguous way to study the deconhnement 
transition; one investigates the chiral transition. Although 
commonly treated to be equivalent, there is no reason why, 
or if at all, these two phase transitions should be simulta¬ 
neous or of the same order (Alam et al. 1996). The order of 
the chiral phase transition depends rather crucially on the 
strange quark mass. If the strange quark is heavy, then the 
chiral phase transition is probably of first order. Otherwise, 
it may be of second order. The strange quark mass being of 
the order of the QCD scale, the situation is still controver¬ 
sial (Blaizot 1999). There are additional ambiguities arising 
from finite size effects of the lattice which may tend to mask 
the true order of the transition. Our interest here is in the 
deconfinement transition. If it is indeed of hrst order, the 
Hnite size effects which could mask it would be negligibly 
small in the early universe. In such circumstances, Witten 
(1984) argued, in a seminal paper, that strange quark matter 
could be the true ground state of Quantum Chromodynam¬ 
ics (QCD) and that a substantial amount of baryon number 
could be trapped in the quark phase which could evolve into 
strange quark nuggets (SQNs) through weak interactions. 
(For a brief review of the formation of SQNs, see Alam, 
Raha & Sinha (1999).) QCD - motivated studies of baryon 
evaporation from SQNs have established (Bhattacharjee et 
al. 1993; Sumiyoshi & Kajino 1991) that primordial SQNs 
with baryon numbers above ~ ]^q 40-42 cosmologi¬ 

cally stable. We have recently shown that without much hne 
tuning, these stable SQNs could provide even the entire clo¬ 
sure density (H ~ 1) (Alam et al. 1999). Thus, the entire 


cold dark matter (CDM) ( Hcdm ~ 0.3-0.35 ) could easily 
be explained by stable SQNs. 

We can estimate the size of the SQNs formed in the first 
order cosmic QCD transition in the manner prescribed by 
Kodama, Sasaki and Sato (Kodama et al. 1982) in the con¬ 
text of the GUT phase transition. For the sake of brevity, 
let us recapitulate very briefly the salient points here; for 
details, please see Alam et al (1999) and Bhattacharyya et 
al (2000). Describing the cosmological scale factor R and 
the coordinate radius X in the Robertson-Walker metric 
through the relation 

ds^ = —dt^ dx^ 

= -dt^-\-R^{dX^-\-X^{sinedcj>'^-\-de^)}, ( 1 ) 


one can solve for the evolution of the scale factor R{t) in the 
mixed phase of the first order transition. In a bubble nucle- 
ation description of the QCD transition, hadronic matter 
starts to appear as individual bubbles in the quark-gluon 
phase. With progressing time, they expand, more and more 
bubbles appear, coalesce and finally, when a critical fraction 
of the total volume is occupied by the hadronic phase, a 
continuous network of hadronic bubbles form (percolation) 
in which the quark bubbles get trapped, eventually evolv¬ 
ing to SQNs. The time at which the trapping of the false 
vacuum (quark phase) happens is the percolation time tp, 
whereas the time when the phase transition starts is de¬ 
noted by ti. Then, the probability that a spherical^^ region 
of co-coordinate radius X lies entirely within the quark bub¬ 
bles would obviously depend on the nucleation rate of the 
bubbles as well as the coordinate radius X{tp, ti) of bubbles 
which nucleated at ti and grew till tp. For a nucleation rate 
I{t), this probability P{X,tp) is given by 


P{X,tp) = exp 


dvr 

'Y 


dtI{t)R\t)[X + X{tp,ti)f 


■ ( 2 ) 


After some algebra (Bhattacharyya et al. 2000), it can be 
shown that if all the cold dark matter (CDM) is believed to 
arise from SQNs, then their size distribution peaks, for rea¬ 
sonable nucleation rates, at baryon number ~ 10'*^“^^, evi¬ 
dently in the stable sector. It was also seen that there were 
almost no SQNs with baryon number exceeding 10^®“^^, 
comfortably lower than the horizon limit of ~ 10®® baryons 
at that time. Since Hb is only about 0.04 from BBN, Hcdm 
in the form of SQNs would correspond to ~ 10®^ baryons so 
that there should be 10^“® such nuggets within the horizon 
limit at the microsecond epoch, just after the QCD phase 
transition (Alam et al. 1998; Alam et al. 1999). We shall 
return to this issue later on. 

It is therefore most relevant to investigate the fate of 
these SQNs. Since the number distribution of the SQNs is 
sharply peaked (Bhattacharyya et al. 2000), we shall assume, 
for our present purpose, that all the SQNs have the same 
baryon number. 

The SQN’s formed during the cosmic QCD phase tran¬ 
sition at T ~ 100 MeV have high masses (~ 10'*'*GeV) and 
sizes {Rn ~ Im) compared to the other particles (like the 
usual baryons or leptons) which inhabit this primeaval uni¬ 
verse. These other particles cannot form structures until the 


tt For the QCD bubbles, it is believed that there is a sizable 
surface tension which would facilitate spherical bubbles. 
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temperature of the ambient universe falls below a certain 
critical temperature characteristic of such particles; till then, 
they remain in thermodynamic equilibrium with the radia¬ 
tion and other species of particles. This characteristic tem¬ 
perature is called the freezeout temperature for the corre¬ 
sponding particle. Obviously the freezeout occurs earlier for 
massive particles for the same interaction strength. In the 
context of cosmological expansion of the universe this has 
important implications ; the ’frozen’ objects can form struc¬ 
tures. These structures do not participate in the expansion 
in the sense that the distance between the subparts do not 
increase with the scale size and only their number increases 
due to the cosmological scale factor. 

For the SQN’s, however, the story is especially inter¬ 
esting. Even if they continue to be in kinetic equilibrium 
due to the radiation pressure (photons and neutrinos) acting 
on them, their velocity would be extremely non relativistic. 
Also their mutual separation would be considerably larger 
than their radii; for example, at T ~ lOOMeV, the mutual 
separation between the SQN’s (of size ~ 10'*'* baryons) is 
estimated to be around ~ 300m. It is then obvious that the 
SQN’s do not lend themselves to be treated in a hydrody- 
namical framework; they behave rather like discrete bodies 
in the background of the radiation fluid. They thus expe¬ 
rience the radiation pressure, quite substantial because of 
their large surface area as well as the gravitational potential 
due to the other SQN’s. 

In such a situation, one might be tempted to assume 
that since the SQN’s are distributed sparsely in space and 
interact only feebly with the other SQN’s through gravita¬ 
tional interaction, they might as well remain forever in that 
state. This, in fact, is quite wrong, as we demonstrate below. 

The fact that the nuggets remain almost static is hardly 
an issue which requires justification. The two kinds of mo¬ 
tion that they can have are random thermal motion and 
the motion in the gravitational well provided by the other 
SQN’s. This other kind of motion is typically estimated us¬ 
ing the virial theorem, treating the SQN’s as a system of par¬ 
ticles moving under mutual gravitational interaction (Bhatia 
2001; Peebles 1980). The kinetic energy { K ) and potential 
energy V of the nuggets at temperature T = 100 MeV can 
be estimated as, 


I/ = ^G 


K = ^NkbT 

MiMj GM^N^ 
Ri,j 2R av 


(3) 


where kb is the Boltzmann constant. Mi, Mj are the masses 
of the ith and jth nugget, Ri,j is the distance between them 
and Rav is the average inter-nugget distance. Substituting 
the number of nuggets N — 10^, the baryon number of each 
nugget to be 10** and Rav = 300m, one gets K = 2.4 x 10“* 
and V — 3.09 x 10^® (in MKS units) so that the ratio of K 
and y becomes ~ 10”®®. Thus it is impossible for these 
objects to form stable systems, orbiting round each other. 
On the other hand the smallness of the kinetic energy shows 
that gravitational collapse might be a possible fate. 

Such, of course, would not be the case for any other mas¬ 
sive particles like baryons; their masses being much smaller 
than SQN, the kinetic energy would continue to be very 
large till very low temperatures. More seriously, the Virial 


theorem can be applied only to systems whose motion is sus¬ 
tained. For SQNs, a notable property is that they become 
more and more bound if they grow in size. Thus SQNs would 
absorb baryons impinging on them and grow in size. Also, 
if two SQNs collide, they would naturally tend to merge. In 
all such cases, they would lose kinetic energy, making the 
Virial theorem inapplicable. 

One can argue that the mutual interaction between uni¬ 
formly dispersed particles would prevent these particles from 
forming a collapsed structure, but that argument holds only 
in a static and infinite universe, which we know our universe 
is not. Also a perfectly uniform distribution of discrete bod¬ 
ies is an unrealistic idealization and there must exist some 
net gravitational attraction on each SQN. The only agent 
that can prevent a collapse under this gravitational pull is 
the radiation pressure, and indeed its effect remains quite 
substantial until the drop in the temperature of the ambi¬ 
ent universe weakens the radiation pressure below a certain 
critical value. In what follows, we try to obtain an estimate 
for the point of time at which this can happen. 

It should be mentioned at this juncture that for the 
system of discrete SQN’s suspended in the radiation fluid, a 
detailed numerical simulation would be essential before any 
definite conclusion about their temporal evolution can be 
arrived at. This is a quite involved problem, especially since 
the number of SQN’s within the event horizon, as also their 
mutual separation, keeps increasing with time. Our purpose 
in the present work is to examine whether such an effort 
would indeed be justified. 

Let us now consider the possibility of two nuggets co¬ 
alescing together under gravity, overcoming the radiation 
pressure. The mean separation of these nuggets and hence 
their gravitational interaction are determined by the tem¬ 
perature of the universe. If the entire CDM comes from 
SQNs, the total baryon number contained in them within 
the horizon at the QCD transition temperature ( ~ 100 
MeV) would be ~ 10®* (see above). For SQNs of baryon 
number fejv each, the number of SQNs within the horizon 
at that time would be just (l0®*/6jv). Now, in the radia¬ 
tion dominated era the temperature dependence of density 
un ~ r®, horizon volume Vh varies with time as t®, i.e. 
Vh ~ T”® and hence the variation of the total number in¬ 
side the horizon volume will be Nn ~ T”®. So at any later 
time, the number of SQNs within the horizon ( Nm ) and 
their density ( un ) as a function of temperature would be 
given by : 


Nn{T) 

nN{T) 


10®* /100MeV\® 
ITT V T ) 
Nn 3Viv 
VT ^ 47r(2t)3 


(4) 

(5) 


where the time t and the temperature T are related in the 
radiation dominated era by the relation : 


t = 0.3(? 


-1/2 Wlpi 
* ^2 


( 6 ) 


with (/* being ~ 17.25 after the QCD transition (Alam et al. 
1999). 

From the above, it is obvious that the density of SQNs 
decreases as t”®/® so that their mutual separation increases 
as t*/®. Therefore, the force of their mutual gravitational pull 
will decrease as t”*. On the other hand, the force due to the 
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radiation pressure (photons and neutrinos) resisting motion 
under gravity would be proportional to the radiation energy 
density, which decreases as or It is thus reasonable 
to expect that at some time, not too distant, the gravita¬ 
tional pull would win over the radiation pressure, causing 
the SQNs to coalesce under their mutual gravitational pull. 
The expression for the gravitational force as a function of 
temperature T can written as : 




(7) 


where fojv is the baryon number of each SQN and m„ is the 
baryon mass. f„„(T) is the mean separation between two 
nuggets and is given by the cube root of the ratio k of total 
volume available and the total number of nuggets 

1.114 X 


The force due to the radiation pressure on the nuggets may 
be roughly estimated as follows. We consider two objects 
(of the size of a typical SQN) approaching each other due to 
gravitational interaction, overcoming the resistance due to 
the radiation pressure. The usual isotropic radiation pres¬ 
sure is |pc^, where p is the total energy density, includ¬ 
ing all relativistic species. The nuggets will have to over¬ 
come an additional pressure resisting their mutual motion, 
which is given by ipc^( 7 — 1 ); the additional pressure arises 
from a compression of the radiation fluid due to the mo¬ 
tion of the SQN. The moving SQN would become a pro¬ 
late ellipsoid (with its minor axis in the direction of motion 
due to Lorentz contraction), whose surface area is given by 
'^ ). The eccentricity e is related to the Lorentz 

factor 7 as e = ^ . For small values of e (small 7 ), 

sin~^e ~ e, so that the surface area becomes . 

Thus the total radiation force resisting the motion of SQNs 
is 

1 2 

74ad = -pradCUfall(7r7?jv)/37 (9) 


where prad is the total energy density at temperature T, Wfaii 
or /3c is the velocity of SQ Ns dete rmined by mutual grav¬ 
itational field and 7 is l/y^l — The quantities Fiad, /3 
and 7 all depend on the temperature of the epoch under 
consideration. (It is worth mentioning at this point that the 
t dependence of Fi-ad is actually sharper than the 

estimated above, because of the Ufaii, which goes as 
The ratio of these two forces is plotted against tempera¬ 
ture in figure 1 for two SQNs with initial baryon number 

each. It is obvious from the figure that ratio F’grav/F/ad 
is very small initially. As a result, the nuggets will remain 
separated due to the radiation pressure. For temperatures 
lower than a critical value Tci, the gravitational force starts 
dominating, facilitating the coalescence of the SQNs under 
mutual gravity. 

Let us now estimate the mass of the clumped SQNs, 
assuming that all of them within the horizon at the critical 
temperature will coalesce together. This is in fact a conser¬ 
vative estimate, since the SQNs, although starting to move 
toward one another at Td, will take a finite time to actually 
coalesce, during which interval more SQNs will arrive within 
the horizon. 

In table 1, we show the values of Td for SQNs of dif- 



Figure 1. Variation of the ratio Tgrav/Trad with temperature. 
The dot represents the point where the ratio assumes the value 
1 . 


Table 1. Critical temperatures ( Td ) of SQNs of different initial 
sizes fejV) the total number Nff of SQNs that coalesce together 
and their total final mass in solar mass units. 


fejv 

Td (MeV) 

Nn 

M/MQ 

10^2 

1.6 

2.44 X 10^"^ 

0.24 

1044 

4.45 

1.13 X IQii 

0.01 

1046 

20.6 

1.1 X 10 ’’ 

0.0001 


ferent initial baryon numbers along with the final masses of 
the clumped SQNs under the conservative assumption men¬ 
tioned above. 

It is obvious that there can be no further clumping of 
these already clumped SQNs; the density of such objects 
would be too small within the horizon for further clump¬ 
ing. Thus these objects would survive till today and perhaps 
manifest themselves as MACHOs. It is to be reiterated that 
the masses of the clumped SQN’s given in table 1 are the 
lower limits and the final masses of these MACHO candi¬ 
dates will be larger. (The case for 6 jv = 10"^® is not of much 
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interest, especially since such high values of bN are unlikely 
for the reasonable nucleation rates (Bhattacharyya et al. 
1999; Bhattacharyya et al. 2000); we therefore restrict our¬ 
selves to the other cases in table 1 in what follows.) A more 
detailed estimate of the masses will require a detailed sim¬ 
ulation, but very preliminary estimates indicate that they 
could be 2-3 times bigger than the values quoted in table 1. 

The total number of such clumped SQNs ( Amacho ) 
within the horizon today is evaluated in the following way. 
With the temperature ~ 'PK and time ~ 4 x 10^^ sec¬ 
onds, the total amount of visible baryons within the hori¬ 
zon volume can be evaluated using photon to baryon ratio 
rj ~ 10“^°. The amount of baryons in the CDM will be 
times the total number of visible baryons. This comes out 
to be ~ 1.6 X 10™, flcDM and fin being 0.3 and 0.01 respec¬ 
tively. The total number of baryons in a MACHO is x Nn 
i.e. 2.44 X 10®® and 1.13 x 10®® for initial nugget sizes 10"^^ 
and 10^^ respectively. The quantities &iv and Ajv are taken 
from the Table 1. So dividing the total number of baryons 
in CDM by that in a MACHO, the Amacho comes out to be 
in the range ~ 10^®“^"'. 

We can also mention here that if the MACHOs are in¬ 
deed made up of quark matter, then they cannot grow to 
arbitrarily large sizes. Within the (phenomenological) Bag 
model picture (Chodos et al. 1974) of QCD confinement, 
where a constant vacuum energy density (called the Bag 
constant) in a cavity containing the quarks serves to keep 
them conhned within the cavity, we have earlier investigated 
(Banerjee et al. 2000) the upper limit on the mass of astro- 
physical compact quark matter objects. It was found that 
for a canonical Bag constant B of (145 MeV) this limit 
comes out to be 1.4 Mq. The collapsed SQNs are safely be¬ 
low this limit. (It should be remarked here that although 
the value of B in the original MIT bag model is taken to 
be B = 145 MeV from the low mass hadronic spectrum, 
there exist other variants of the Bag model (Hasenfratz & 
Kuti 1978) , where higher values of B are required. Even 
for B = 245 MeV, this limit comes down to 0.54 Mq 
(Banerjee et al. 2000), which would still admit such SQN. 

As a consistency check, we can perform a theoretical 
estimate of the abundance of such MACHO’s in the galactic 
halo which is conventionally given by the optical depth. The 
optical depth is the probability that at any instant of time a 
given star is within an angle 9e of a lens, the lens being the 
massive body (in our case MACHO) which causes the de¬ 
flection of light. In other words, optical depth is the integral 
over the number density of lenses times the area enclosed 
by the Einstein ring of each lens. The expression for optical 
depth can be written as (Narayan & Bartelmann 1999): 

47rG 9 f 

T = ——Dg / p{x)x{l — x)dx (10) 

where Ds is the distance between the observer and the 
source, G is the gravitational constant and x = DdDs~^, 
Dd being the distance between the observer and the lens. In 
particular p is the mass-density of the MACHOs, which is of 
the form p ~ po^ in the naive spherical halo model, which 
we have adopted in our calculations. In the present case po 
is given by 

Afmacho X N macho /1 1 \ 


where R — + Dg + 2DeDs cos<(), (j) and Dg being the 

inclination of the LMC and the distance of observer (earth) 
from the Galactic centre respectively. Mmacho and Amacho 
are the mass of a MACHO and the total number of MACHOs 
in the Milky Way halo. 

The total visible mass of the Milky Way ( ~ 1.6 x 
10®® Mq ) corresponds to ~ 2 x 10®® baryons. This corre¬ 
sponds to a factor of ~ 2 x 10“® of all the visible baryons 
within the present horizon. Scaling the number of clumped 
SQNs within the horizon by the same factor yields a to¬ 
tal number of MACHOs, Amacho ~ 10®®“®® in the Milky 
Way halo for the range of baryon number of initial nuggets 
&JV = 10®^“®®. The value of Dg and Dg are taken to be 10 
and 50 kpc, respectively. The value of the inclination angle 
used here is 40 degrees. Using these values for a naive in¬ 
verse square spherical model comprising such objects upto 
the LMC, we obtain an optical depth of ~ 10“® — 10“®". The 
uncertainty in this value is mainly governed by the value 
of 77 , Hcdm and Hb, and to a lesser extent by the specific 
halo model. This value compares reasonably well with the 
observed value and may be taken as a measure of reliability 
in the proposed model. 

As an interesting corollary, let us mention that the sce¬ 
nario presented here could have other important astrophys- 
ical signihcance. The origin of cosmic rays of ultra-high en¬ 
ergy > 10®® eV continues to be a puzzle. One of the pro¬ 
posed mechanisms (Bhattacharjee & Sigl 2001) envisages a 
top-down scenario which does not require an acceleration 
mechanism and could indeed originate within our galactic 
halo. For our picture, such situations could easily arise from 
the merger of two or more such MACHOs, which would shed 
the extra matter so as to remain within the upper mass limit 
mentioned above. This is currently under active investiga¬ 
tion. 

We thus conclude that gravitational clumping of the pri¬ 
mordial SQNs formed in a first order cosmic quark - hadron 
phase transition appears to be a plausible and natural expla¬ 
nation for the observed halo MACHOs. It is quite remark¬ 
able that we obtain quantitative agreement with the exper¬ 
imental values without having to introduce any adjustable 
parameters or any hue-tuning whatsoever. We may hnish 
by quoting a famous teaching of John Archibald Wheeler, 
“One should never do a calculation unless one knows an an¬ 
swer” . Our attempt in this work has been to hnd an answer 
so that a calculation (in this case, a detailed simulation of 
the collapsing SQNs) can be embarked upon. 
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their warm hospitality. 


REFERENCES 

Alam J., Raha S., Sinha B., 1996, Phys. Rep., 273, 243 

Alam J., Raha S., Sinha B., 1998, Nucl. Phys. A638, 523c 

Alam J., Raha S., Sinha B., 1999, ApJ, 513, 572 

Alcock C. et ah, 1993, Nat, 365, 621 

Alcock C. et ah, 2000, ApJ, 542, 281 

Aubourg E. et ah, 1993, Nat, 365, 623 

Aubourg E. et ah, 1999, A&A, 347, 850 

Banerjee S., Ghosh S. K., Raha S., 2000, J. Phys., G26, LI 


© 0000 RAS, MNRAS 000, 000-000 



6 Banerjee et al. 


Bhatia V. B., 2001, A textbook on Astronomy and Astrophysics, 
with elements of Cosmology, Narosa Publishing, New Delhi 
Bhattacharjee P., Alam J., Sinha B., Raha S., 1993, Phys. Rev., 
D48, 4630 

Bhattacharjee P., Sigl G., 2000, Phys. Rep., 327, 109 
Bhattacharyya A., Alam J., Sarkar S., Roy P., Sinha B., Raha S., 
Bhattacharjee P., 1999, Nucl. Phys., A661, 629c 
Bhattacharyya A., Alam J., Sarkar S., Roy P., Sinha B., Raha S., 
Bhattacharjee P., 2000, Phys. Rev., D61, 083509 
Blaizot J. P., 1999, Nucl. Phys., A661, 3c 

Chodos A., Jaffe R. L., Johnson K., Thorn C.B., Weisskopf V. 

F., 1974, Phys. Rev., D9, 3471 
Evans N. W., Gyuk G., Turner M. S., 1998, ApJ, 501, L45 
Fields B. D., Freese K, Graff D. S., 1998, New Astron., 3, 347 
Freese K., Fields B. D., Graff D. S., 2000, astro-ph/0002058 
Gyuk G., Gates E., 1999, MNRAS, 304, 281 
Hasenfratz P., Kuti J., 1978, Phys. Rep., 40, 75 
Jedamzik K., 1998, Phys. Rep., 307, 155 
Jedamzik K., Niemeyer J.G., 1999, Phys. Rev., D59, 124014 
Kodama H., Sasaki M. and Sato K., 1982, Prog. Theo. Phys., 8, 
1979 

Narayan R., Bartelmann M., 1999, in Dekel A., Ostriker J. P., eds. 
Formation of structure in the universe. Cambridge University 
Press, Cambridge 
Paczynski B., 1986, ApJ, 304, 1 

Peebles P. J. E., 1980, The large scale structure of the universe. 

Princeton Univ. Press, Princeton, NJ 
Schmid C., Schwarz D. J., Widerin P., 1999, Phys. Rev., D59, 
043517 

Schramm D. N.,1998, in Sinha B., Viyogi Y. P., Srivastava D. K., 
eds, Proc. ICPA-QGP’97, Physics and Astrophysics of Quark- 
Gluon Plasma. Narosa Publishing, New Delhi, p. 29 
Sumiyoshi K., Kajino T., 1991, Nucl. Phys. (Proc. Suppl.), 24, 80 
Sutherland W, 1999, Rev. Mod. Phys., 71, 421 
Turner M. S., 1999, Rev. Mod. Phys., 71, S145 
Turner M. S., 2000, Phys. Scr., T85, 210 
Witten E., 1984, Phys. Rev., D30, 272 
Zaritsky D., Lin D. N. C., 1997, AJ, 114, 2545 


© 0000 RAS, MNRAS 000, 000-000 



